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A semiclassical method is used to study Landau damping of transverse pseudo- 
spin waves in harmonically trapped ultracold gases in the collisionless Boltz- 
mann limit. In this approach, the time evolution of a spin is calculated nu- 
merically as it travels in a classical orbit through a spatially dependent mean 
field. This method reproduces the Landau damping results for spin-waves in 
unbounded systems obtained with a dielectric formalism. In trapped systems, 
the simulations indicate that Landau damping occurs for a given spin-wave 
mode because of resonant phase space trajectories in which spins are "kicked 
out" of the mode (in spin space). A perturbative analysis of the resonant and 
nearly resonant trajectories gives the Landau damping rate, which is calcu- 
lated for the dipole and quadrupole modes as a function of the interaction 
strength. The results are compared to a numerical solution of the kinetic 
equation by Nikuni et al. 



1. INTRODUCTION 

Recent experiments at JILA 1 2 have demonstrated transverse spin waves 
in the pseudo spin dynamics of the hyperfine levels of trapped ultra-cold 
87 Rb atoms in the nondegenerate regimes. As in spin-polarized Boltzmann 
gases, an exchange mean-field supporting collective spin behavior can occur 
when the thermal de Broglie wavelength exceeds the effective range of the 
interaction between colliding atomP^. Numerical simulations of the appro- 
priate kinetic equations with JILA parameters have been shown to reproduce 
the observed pseudo-spin oscillations.^ Using a moments method, Nikuni et 
a/P obtained analytic expressions for the spin-wave frequencies and damp- 
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ing rates of the dipole and quadrupole modes in a trapped Boltzmann gas. 
The results were found to be in excellent agreement with a one dimensional 
numerical solution of the kinetic equation, except at intermediate densities, 
where Landau damping was thought to be important. However, an estimate 
using the Landau damping result for a homogeneous syste m 7 yielded a rate 
that was order of magnitude too large. 

One of the disadvantages of the moments method is that the truncated 
form of the distribution functions do not account for the mean field coupling 
of the collective modes to higher excitations that give rise to Landau damp- 
ing. Indeed, Landau damping usually requires a pole in the spin distribution 
functionP A trapped system poses another problem in that the excitation 
spectrum is discrete, whereas Landau damping calculations usually involve 
integration over a continuous spectrum. One method that is often used is to 
look at elongated clouds and treat the long axis of the trap as infinite p but 
the applicability of this approximation to low lying states is questionable. 
In the following we use a semi-classical approach to extend the moments 
method to account for Landau damping in trapped systems. 



2. HOMOGENEOUS SYSTEMS 

Before considering trapped systems, we review Landau damping in a 
homogeneous system. We start with the kinetic equation for the normalized 
transverse spin distribution s + (z,p) in the collisionless regime, linearized for 
small tip angles 

d t s + (z,p) = -(p/m)d z s + (p,z) -iQ[s + (z,p) - S + (z)] (1) 

where ft is the mean-field frequency, S + (z) = J s + (z , p)n(p) dp and n(p) = 
(27rm/cBT)~ 1//2 exp[— p 2 /(2m/csT)] is the normalized distribution in momen- 
tum space. Now substituting traveling wave solutions of the form s + (z,p) = 
s + (p) exp[i(qz — ut)], and solving for s + (p) we get 

s +(p) = QS+ (2) 

K ; qp/m + 

Integrating over p and taking into account the pole at pl = (m/q)(u> — ft) 
we get the same result obtained in Ref.7 using an RPA approximation: 

ksT 2 ft 2 I nm 



in 



2 q 2 k B T 



(3) 
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3. SEMICLASSICAL APPROACH 

The above calculation of the damping rate is direct and simple, but it 
cannot be easily adapted to trapped systems. We now recalculate the above 
damping rate for the unbounded case using a semiclassical approach. In 
the next section we adapt this approach to a numerical method for trapped 
systems. 

If we consider the small transverse component of a spin as it travels 
along the helix, we find that its equation of motion is given by da(t)/dt = 
Qa(t) x S[z(t)]. In linearized form this becomes 

d + = -in (a + - S + [z{t)}) = -m (a + - e nz{t) e -i"A (4) 

where in the collisionless regime, u> q = —kBT/Qmq 2 and z(t) = pt/m. The 
solution is given by 

a + (t) = 2 ,n Sin[(gp/m + (1 - u q )t/2} ci(qp /m-n-u a )t/2 (5) 
V ; (qp/m + Q-lu) ' V ; 

Obviously, we have resonance at p^ = (tu q — £l)m/q in which case we have 
= Qt. Otherwise, a + is oscillatory, with an amplitude 2mQ/q(p—pL) 
and a beat frequency oJbeat = q(p — Pl)/ 171 - We note that in the absence of 
damping the amplitude of a + (j~ is constant. Thus, the derivative of <t + £t~ 
gives the Landau damping rate. Using Eq|5]we have 

lim \a{t)\ 2 /t = 2-ir 1 ^-5(p - {uj q ~ tym/q) (6) 
t — ^00 q 

The contribution to the Landau damping of s + is half this valueP To find 
the total landau damping rate we integrate over the p distribution according 
to 

7L = 7r n(p)5(p - p L ) dp (7) 



which gives, assuming w«Sl, the same result as EqO 

We now outline a numerical procedure to find the Landau damping 
rate that can be generalized to trapped systems. First we search for the 
resonant momentum p^. We then measure the slope of the growth of |<7 + (£)| 
at resonance to get Q e ff- Finally we find q e ff from the beat frequency 
exhibited by s + (t) when p is slightly off resonance. According to EqEJwc 
have q e ff(p~ PL)T beat /m = 2ir or 
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4. TRAPPED SYSTEMS 

We restrict our attention to elongated (quasi-ld) clouds. For a trapped 
system the trajectories through phase space are no longer straight lines of 
constant p, but ellipses of constant energy E = p 2 /2m + muj 2 z 2 /2, where lo z 
is the axial trap frequency. Eq|I] generalizes according to 



a 



-iQ(z)(a + - S + {z,t)) (9) 



with z(t) = \J2E/ (mu 2 ) cos(co z t) and Q(z) = Qon(z)/no, where Q,q is ra- 
dially averaged exchange frequency^ and n(z)/no = exp(— muj 2 z 2 /2ksT) is 
the ratio of the density of the cloud to its peak value. For the transverse 
magnetization we use as an approximation the results from the moments 
method of Nikuni et for the dipole and quadrupole modes: 

S+(z) = Cl ze- iuit (10) 
S+(z) = c 2 {z 2 -l)e~^ (11) 

where the frequencies are given by uij = -n /y/E + y/SF/8 + j 2 u;'i and the 
normalization constants Cj are determined from the condition 

1 = J \s + (z,p)\ 2 n(z,p)dzdp (12) 

where n(z,p) = uj z /(2irkBT)exp(—E(z,p)/kT). The Landau damping rate 
is then given by 

1l = ~ [ ^—j^-n(z,p)dzdp 

2E ' e- E ' kBT 5{E-E L )dE 




(13) 

To calculate the Landau damping rate we use the numerical proce- 
dure described in the previous section. Solutions of EqEl (See Fig.l) show 
the same sort of resonant behavior as the exact solution for the homoge- 
neous system. However, the effective parameters for the trapped system are 
very different from those of the homogeneous system. For f^o = 2w 2 in the 
quadrupole mode we find a much smaller effective value Q e ff = 0.38o; z . For 
the same fio we get q e ff = Q^lz^ 1 where zr = (ksT /muj 2 ) 1 / 2 is the ther- 
mal cloud size. This is much smaller than an estimate in which the distance 
between the nodes is taken to be a half a wavelength, giving q = 1.57lo z . (Not 
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Fig. 1. Time evolution of the transverse component of a spin as it moves 
through the dipole mode for = 3.5uj z at: the resonant energy E = El 
(straight line), slightly off resonance showing beats, and far from resonance. 



surprisingly, the dipole mode has a smaller value q e ff = 0.158z^ .) Again 
for Qq = 2uj z , in the trapped system we find El = 4.36/cbT whereas the 
argument of the exponential in Eq|3] gives 0.811. Fig. 2a shows the Landau 
damping rate calculated from Eq ll.SI Fig. 2b shows the collisional damping 
rate calculated via the moments methocP as well as the sum of the Landau 
and collisional damping rates. For comparison the numerical solutions of 
the kinetic equation of Nikuni et are also shown. Although the location 
of the Landau damping peaks are correct and they are of the right order 
of magnitude, there are qualitative differences with the numerical solutions. 
The peak Landau damping rate for the dipole mode is too large, whereas it 
is too large for the quadrupole at high densities. 

There are two possible improvements to the model that immediately 
come to mind: 1) Better approximations of the functions S^(z) and the 
frequencies ojj could be used. These could be obtained by including higher 
moments in the moments calculation and should give a more accurate value 
of the exchange field experienced by spins, especially in orbits with E > ksT. 
2) The effects of collisions on Landau damping could be included. Although 
these might be rather difficult to incorporate within the framework of the 
current model, they might be easy to analyze for the homogeneous system, 
in Eqs.(l-3), for example, and perhaps the predictions of the trapped model 
could then be modified in a simple way. We leave this for future work. 
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Fig. 2. a) Landau damping rate of the dipole (dashed) and quadrupole (solid) 
modes as a function of the exchange frequency, b) Damping rates calculated 
by Nikuni et with a numerical solution of the kinetic equation for the 
dipole (circles) and quadrupole (squares) modes for 87 Rb at T =800nK. 
Also shown are the damping rates calculated with the moments methocP 
(dashed), and with Landau damping included (solid). 
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